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ABSTRACT 

A new method is presented for the construction of a nat- 
ural continuous wavelet transform on the sphere. It incorpo- 
rates the analysis and synthesis with the same wavelet and the 
definition of translations and dilations on the sphere through 
the spherical harmonic coefficients. We construct a couple of 
wavelets as an extension of the flat Mexican Hat Wavelet to 
the sphere and we apply them to the detection of sources on 
the sphere. We remark that no projections are used with this 
methodology. 

1. INTRODUCTION 

Multiscaling analysis techniques dealing with the analy- 
sis/synthesis of nD-images defined on intervals of have 
been applied in many fields of physics in the last 15 years. 

For instance, in the case n = 1 one has electronics and 
audio signals, in the case n = 2 one has optical or infrared 
images whereas for n = 3 one deals with fluid dynamics or 
the large-scale structure of the universe as 3D-images. How- 
ever, there are data given on other manifolds like the circle 
S\ (e. g. scanning along circles the microwave sky) and the 
sphere (e. g. geophysics). In this paper, we are interested 
in data distributed on the sphere. Trivially, for the study of 
local properties (e. g. detection of objects) one can project on 
the tangent plane at any point on the sphere to make this type 
of analysis but when global properties are taken into account 
the curvature of the sphere can not be neglected. 

A first approach to deal with these global properties is to 
make some global projection of all the points of the sphere. 
The stereographic projection has been recently used deal- 
ing with the continuous wavelet transform. In this case, to 
get the wavelet coefficient at any point on the sphere, one 
projects from the opposite point to the local tangent plane. 
[1] have made a connection to group theory. The transla- 
tions and dilations in the wavelet have their definition on the 
plane. Clearly, such a projection does not take into account 
the topological structure of the sphere. Some applications 
to cosmology, in particular the study of anisotropics of the 
cosmic microwave background radiation have been done by 
some authors ([3], [7], [10]) using the projection of the Mexi- 
can hat wavelet. A drawback of such projection is the obvi- 
ous deformation of the pixels and wavelets near the projec- 
tion pole. We remark that the synthesis can be done in terms 
of another biorthogonal wavelet [11]. 

Another approach uses some analyzing wavelet functions 
that are defined in terms of spherical harmonics [5] with 
a definition of the dilation operator and conditions on the 
wavelets in such a way to get a synthesis formula. The draw- 



back of such methodology is: the dilations do not satisfy 
the appropriate flat limit in general. Also some examples 
of wavelet functions are poorly localized (e. g. Abel-Poison 
wavelets). 

A different approach assumes from the beginning dis- 
crete wavelets incorporating tensor product approaches in 
polar coordinates, then the two poles are singular points re- 
garding approximation/stability properties ([4], [6]). Another 
approach is adapted to arbitrary point systems or triangula- 
tions on the spheres, then there is no efficient tool as fast 
wavelet algorithms. In the approach by [8] basis are defined 
on a quasi-uniform icosahedral triangulation on the sphere al- 
lowing for a fast algorithm. However, biorthogonal wavelets 
are needed and a lifting scheme for the multiresolution is ap- 
plied avoiding the concepts of translations and dilations and 
also it is not clear whether the construction leads to a stable 
^2(52) basis. Haar-type wavelets have been developed us- 
ing different pixel combinations ([2], [7], [9]). The first case 
uses the lifting scheme weighting for the area of the pixels 
whereas in the other two cases an equal area pixelization is 
used but the Haar-type transform is only applied on regions 
of the sphere covering only -j^ of the total area. Clearly, with 
any pixelization the symmetry on the sphere is lost. 

In this paper we will consider a continuous approach, we 
will introduce a methodology that incorporates the analysis 
and synthesis of any function defined on the sphere ^2 us- 
ing the same circularly-symmetric wavelet and also we will 
introduce the generalization of the translations/dilations. In 
this sense we follow Freeden's approach working with spher- 
ical harmonics. Examples will be given that have the appro- 
priate flat limit. Finally, the application to the detection of a 
spot is given, studying the concentration of the wavelet coef- 
ficients. 

2. PROPERTIES OF THE WAVELET 

We will consider a circularly-symmetric filter defined on the 
sphere ^2 



^(H-r;^), 



(1) 



where n is a fixed direction. 7 is another fixed, but arbi- 
trary direction, therefore n ■ 7 will represent a rotation on the 
sphere with respect to the direction n defined by the angle 9 
(cos(0) =n-y). R > Q will represent a dilation, which will 
be defined later on through the spherical harmonics. 
We assume the following properties of the filter: 
(i) the analysis of any function /(«) will be done with the 
wavelets ^{n-f,R), 



(ii) the synthesis of any function /(n) will be done with 
the wavelets coefficients and the wavelets ^(n ■ 7;/?), 

(iii) it will incorporate the definition of translation and 
dilation on the sphere. 

We remark that no assumption about compensation of the 
filter (i. e. / dQ.{n)^{n ■y;R)= 0) and projection from to 
^2 is imposed. 

3. ANALYSIS WITH THE FILTER ^ 

We define the wavelet coefficients associated to the transla- 
tion 7 and dilation R for the function f{n) defined on ^2 



(2) 



w(R,y)= da{n)f{nyV{n-f,R). 



Let us assume the standard decomposition of /(n) in spheri- 
cal harmonics Yi„,{n) 

/(«) = Y.flMn), fl,n = [ da{n)f{n)Y;,„{n). (3) 

Im 

By introducing Eq.(|3|l into Eq.(|3 and taking into account 
that Yim{n) is an orthonormal base of ^2, we obtain 



Alt 



^(^' y) = L(^)/'™^/(^)i'/-(7), (4) 

Im 

where the Legendre coefficients associated to the circularly- 
symmetric filter ^ are given by 



^>{n-f,R) = Y.^>i{R)P,[n-y), 



1, 



dyPi{y)^>[y;R) 
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4. SYNTHESIS WITH THE FILTER ^ 



(5) 



Now, let us show that in order to have a reconstruction equa- 
tion, i. e. /(«) as a functional integral of the wavelet coeffi- 
cients and the wavelet base one can impose the condition 



,2/+l 
' 4n 



(6) 



i. e. "P/ {R) depends on the product IR and satisfies the 
admissibility condition 



f°° dl J, , 



(7) 



where / runs in the interval [0,°°). We remark that the analo- 
gous condition to have a reconstruction on the plane by sub- 
stituting I q, q being the wave number in Fourier space. 
Therefore, the filter ^ -given by Eq. (|5jl- can be rewritten as 

^'(n-f,R)=Y,Wi{R)Pi{n-f)^Y.V{lR)Y,l{n)YU7)- (8) 

/ Im 

Firstly, we remark that the previous equation defines a dila- 
tion on the sphere in terms of dilation of the number / and 
a translation on the sphere in terms of a rotation through the 
spherical harmonics Yi,„{f). We think that such a definition is 
the most natural on the sphere and generalizes the one asso- 
ciated to dilations and translations in the plane R^ via Fourier 
space. 



Secondly, we can write the following equation 
dR 



R 



dQ.{f)w{R,y))^'{n-f,R) 
dR, 4n 2 



(9) 



where we have taken the harmonic expansions for w{R,y) 
and ^(« • y;R). If one wants to have this equation propor- 
tional to Y.iinfimYim{n) — f{n), i.e. to be able to reconstruct 
/(n), then it is obvious that necessarily 



dR , 47r 



— ( ) "VjiR^^C^, 



(10) 



where Cif, 7^ must be a constant. Aparticular solution to 
the previous equation is given by Eq. |6|and the admissibility 
condition. In this case, the synthesis equation can be written 

as 

1 fdR 



m ^c^J TJ dmMR.7)'i'in-f,R). (11) 
and the Equation Q can be rewritten as 



AR,f) =Y^fl,nWmYlm{f)- (12) 

Im 

These equations are the analysis/synthesis counterparts on ^2 
of the corresponding ones on R^. 

5. PROPERTIES OF THE FILTER * 

Let us focus on some of the properties of the filter: 

^ is a compensated filter 

Taking into account Eq. 0, C^f^ <°° implies that 
Z*^ , e > as / ^ 0, i. e. !//(/ = 0) = 0. Now, taking into 
account Eq. (|8|l, one obtains 



/ 



dQ.{n)'¥{n-f,R) = ^(0) = 0, 



(13) 



the filter is compensated (hereinafter wavelet). 
Energy of the wavelet 

Taking into account Eq. ^ and the orthonormal property 
of the spherical harmonics, one obtains 

Jd£l{n)^\n-y-R)=Y,^w\lR). (14) 



Energy of the function f{n) 

Taking into account the standard properties of the spher- 
ical harmonics 



/ / dain)fin) = Y,fL fL ^ flnrfL (15) 



We can also prove the following equivalence 



1 fdR 



C\ft J R 



dQ.{n)w-{R,n). (16) 



6. AN EXAMPLE 



7. CONCLUSIONS 



As an example of the previous ideas we will consider the gen- 
eralization to the sphere of the Mexican Hat wavelet (MHW). 
We focus on the MHW because it is a widely used tool in 
Astronomy, well suited for the detection of pointlike ob- 
jects such as extragalactic sources ([3], [7], [10]), but the same 
ideas can be applied to other wavelet families as well. Two 
natural generalizations of the mother Mexican Hat wavelet 
on the plane are possible 

V/i(0-/V^'', xi/2{l)ocl(l + l)e-^'i'+'\ (17) 

where we have taken a unit width to define the mother 
wavelet. For / ^ 1 both functions approach the MHW 
on the plane. We have represented the harmonic coeffi- 
cients \jfi{lR), V/2('^) for different values of /? = 0.2 x l-ij = 
—2,-1,0, 1,2, as well as the profile of the wavelets on real 
space for the same cases, in Figure ^ The differences be- 
tween wavelets Xj/i and v/2 is shown in Figure|2l 

Next, in order to study the concentration in wavelet space 
we will consider a spot with spherical symmetry placed on 
the north pole, i. e. it is defined by a function f{6). In this 
case, we get for the wavelet coefficients 

wiR,e) = Y.MR)Piicose), w,iR)=fiXi/{lR), 
I 

fi = {i + h rdesmePi{cose)f{o). (i8) 

2 Jo 

Let us now consider a very simple spot defined by a top hat 
f{9)~ l[o.eo]- ^ simple way to test how the wavelets allow 
us to concentrate the information in a few number of coeffi- 
cients is to measure somehow the width of the curve w{R, 9) 
of wavelet coefficients for the spot. An intuitive way to do 
this is to define the "energy" as the integral under the squared 
curve w^{R,9) and to see which is the radius 9e that contains 
a given fraction of the total energy. The smaller 9e is, the 
more concentrated the coefficients are. 

We have performed numerical simulations with a sim- 
ple toy model to see which of the two generalizations of the 
MHW, or concentrates more the coefficients. We have 
placed a top hat spot of size 9o ~ 0.2 rad in the North Pole 
and we have filtered it with the wavelets Yi and using dif- 
ferent scales ranging from /? = 0.1to/? = 1.6 rad. The results 
are shown in the upper panel of Figure|3l We have ploted the 
radius 9e{R) such that 68% of the energy is inside the circle 
of radius 9e{R) as a function of the dilation scale R. As can 
be seen, the wavelet i/i concentrates more the coefficients, 
that is, produces smaller values of 9e- 

Now we repeat the same process but using a Gaussian 
spot instead of a top hat. The Gaussian spot is given by 
f{9) — exp(— 0^/20Q ). This case is interesting since most 
of the detectors that operate in microwave Astronomy exper- 
iments have approximately Gaussian response. The lower 
panel of Figure|3]show the results, that are very similar to the 
top hat case. Again, Xj/i concentrates more the coefficients. 
We have tested the resuts for different levels of concentra- 
tion and we have found that up to 75% the wavelet i/i con- 
centrates more the coefficients. Above 75% this is still true 
for large values of the dilation R, while for small R 1/2 con- 
centrates more the coefficients. This is due to the fact that 
y/i produces more small oscilations in the tail of the curve 
w{R,9) for small 



We have developed a constructive wavelet approach on the 
sphere without any projection from the plane. It is a con- 
tinuous transform that allows the analysis and synthesis of 
any function defined on the sphere and incorporates the con- 
cepts of translation and dilation as generalizations of the ele- 
mentary ones defined on the plane. It is a compensated filter 
that conserves the energy of any function. We have consid- 
ered some natural generalizations of the plane Mexican hat 
wavelet and we have applied them to the detection of a big 
spot. The conclusion is that one of the wavelets (y/i) concen- 
trates more the information than the other one. 
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Figure 1: Harmonic coefficients of the wavelets (top left) and i//2 (top right). In the bottom panels, the dilated wavelets Xj/i (left) and \f/2 
(right) are shown on real space. For all the cases, R = 0.2 x 2-^ 





Figure 2; Difference Xjfi — 1/2 on harmonic (top) and real space (bottom) for the same cases as in FigureQ 




Figure 3: Width of tiie wavelet coefficient curves w{R, 6) for the wavelets Yl snd i|/2 applied to a simple top hat spot of size 9q = 0.2 rad 
placed on the north pole and different values of R (top panel) and for a Gaussian spot of size 6q = 0.2 rad placed on the north pole and 
different values of R (lower panel). The width of the curve is defined as the radius 0^ that contains 68% of the energy of the curve. 



